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Abstract 

In this paper we investigate the nature of the transition from Abelian to non- 
Abelian confinement (i.e. crossover vs. phase transition). To this end we 
consider the basic M = 2 model where non- Abelian flux tubes (strings) were 
first found: supersymmetric QCD with the U(iV) gauge group and Nf = N 
flavors of fundamental matter (quarks). The Fayet-Iliopoulos term £ triggers 
the squark condensation and leads to the formation of non- Abelian strings. 
There are two adjustable parameters in this model: £ and the quark mass 
difference Am. We obtain the phase diagram on the (£, Am) plane. At 
large £ and small Am the worldsheet dynamics of the string orientational 
moduli is described by M = 2 two-dimensional CP(iV — 1) model. We 
show that as we reduce £ the theory exhibits a crossover to the Abelian 
(Seiberg-Witten) regime. Instead of iV 2 degrees of freedom of non-Abelian 
theory now only N degrees of freedom survive in the low-energy spectrum. 
Dyons with certain quantum numbers condense leading to the formation of 
the Abelian Zn strings whose fluxes are fixed inside the Cartan subalgebra 
of the gauge group. As we increase N this crossover becomes exceedingly 
sharper becoming a genuine phase transition at N = oo. 



1 Introduction 



Transition from Abelian to non-Abelian confinement emerged as a central 
question in the current explorations of Yang-Mills theories. By Abelian 
confinement we mean that only those gauge bosons that lie in the Cartan 
subalgebra are dynamically important in the infrared, i.e. at distances of the 
order of the inverse flux tube (string) size. By non-Abelian confinement we 
mean such dynamical regime in which at distances of the flux tube formation 
all gauge bosons are equally important. In supersymmetric Af = 2 Yang- 
Mills theories slightly deformed by a /iTr$ 2 term linear confinement was 
discovered [1], explained by the dual Meissner effect. In the limit of small fi 
amenable to analytic studies [1] confinement is Abelian. It is believed that 
as fi gets large, fj, > A, a smooth transition to non-Abelian confinement takes 
place in the Seiberg-Witten model. 

In non-supersymmetric theories a similar purpose construction, with an 
adjustable parameter, was engineered in [2] (see also [3]) where Yang-Mills 
theories on R 3 x Si were considered. The radius of the compact dimen- 
sion r(Si) was treated as a free parameter. At small r(Si), after a center- 
symmetric stabilization, linear Abelian confinement sets in by virtue of the 
Polyakov mechanism [4]. Then it was argued that the transition from the 
small-r(S'i) Abelian confinement regime to the decompactification limit of 
large r, r(Si) ^> A -1 , where confinement is non-Abelian, is smooth. No 
obvious order parameter that could discontinuously change in passing from 
small to large r(Si) was detected. 

This paper presents our new results on this issue. The nature of the 
transition from the Abelian to non-Abelian regime - i.e. phase transition vs. 
crossover - appears to be non universal. If there is a discrete symmetry on the 
string world sheet and the mode of realization of this symmetry changes in 
passing from the Abelian to non-Abelian regime then these two domains are 
separated by a phase transition. A particular nonsupersymmetric example 
[5, 6] of such a situation will be discussed in the bulk of the paper. 

On the other hand, if the mode of realization of the discrete symmetry 
does not change, or there is no appropriate symmetry whatsoever, then the 
Abelian confinement is separated from the non-Abelian one by a crossover 
rather than phase transition. In this paper we focus on Af = 2 Yang-Mills 
theory with the gauge group U(A r ) and Nf = N flavors. To simplify our 
discussion we mostly consider the N = 2 case of U(2) theory with two fla- 
vors. We will show that this model belongs to the second class, with smooth 
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transitions of the crossover type. Although there is no phase transition, we 
find that both perturbative and non-perturbative low energy spectra of the 
theory are drastically changed when we pass from the Abelian to non-Abelian 
regime. 

The benchmark model we will deal with - the U(2) theory with Nf — 2 
(s) quark multiplets - is described in detail in the review paper [7]. It is worth 
recalling that the model is characterized by two adjustable parameters: the 
coefficient of the Fayet-Iliopoulos (FI) term £ and the difference Am of the 
mass terms of the first and second flavors. If £ 3> A 2 , where A is the dy- 
namical scale of the gauge theory at hand, the theory is at weak coupling 
and can be exhaustively analyzed using quasiclassical methods. The domain 
of large \Am\ is that of the Abelian confinement. At small \Am\ confine- 
ment is non-Abelian. A discrete Z2 symmetry inherent to the Lagrangian 
of the world-sheet theory is spontaneously broken in both limits, albeit the 
order parameters are different. Thus we expect (and, in fact, demonstrate) 
a crossover in \Am\. 

The domain of small £ was not considered previously in the context of 
the problem we pose. Our task is to include it in consideration. Various 
regimes of the theory in the (£, Am) plane are schematically shown in Fig. 1. 
We choose Am real, which is always possible to achieve by an appropriate 
U(l) rotation. The vertical axis in this figure denotes the values of the FI 
parameter £ while the horizontal axis represents the quark mass difference. 

As was mentioned, the domain I is that of non-Abelian confinement. 
In this domain the perturbative spectrum of the bulk theory has N 2 light 
states. 1 In the limit of degenerate quark masses the bulk theory has an 
unbroken global SU(2)c + .p symmetry (the so-called color-flavor locking, see 
Sect. 2), and the light states come in adjoint and singlet representations 
of this group. The nonperturbative spectrum contains mesons built from 
monopole-antimonopole pairs connected by two strings, [7]. These strings 
are non-Abelian [8, 9, 10, 11] (see also the reviews [12, 7, 13, 14]). The non- 
Abelian SU(2) part of their fluxes is determined by moduli parameters, whose 
dynamics is described by M = 2 supersymmetric CP(1) on the string world 
sheet. Due to large quantum fluctuations in the CP(1) model the average 
non-Abelian flux of such a string vanishes. 

The domain II is that of Abelian confinement at weak coupling. As 

1 By light we mean those states whose masses are less than or of the order of the inverse 
size of the string. 
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Figure 1: Various regimes in the benchmark J\f 
crossovers. 



= 2 model are separated by 



we increase Am, the W bosons and their superpartners become heavy and 
decouple from the low-energy spectrum. We are left with two photon states 
and their quark M = 2 superpartners. Strings also become Abelian Z 2 
strings. The moduli of the CP(1) model at large Am are fixed in two definite 
directions - the north and south poles of the S 2 sphere (the target space of 
the CP(1) sigma model is S2). The non- Abelian parts of their fluxes no 
longer vanish. 

As we reduce £ and |Am| we enter the domain III. It is nothing but the 
Abelian Seiberg-Witten confinement [1, 15]. The W bosons and their su- 
perpartners decay on the curves of the marginal stability (CMS). The set of 
the light states includes photons and dyons with certain quantum numbers 
(the quarks we started with become dyons due to monodromies as we reduce 
|Am|). Condensation of dyons leads to formation of Abelian Z 2 strings. Non- 
perturbative spectrum still contains mesons built of monopole-antimonopole 
pairs connected by two distinct Z 2 strings. However, now these strings are 
different from those in the domain I. They have nonvanishing fluxes directed 
in the Cartan subalgebra of SU(2). 

These three regimes are separated by crossover transitions. Thus, non- 
Abelian strings can smoothly evolve into Abelian ones and vice versa. This 
is our main finding. 

If, instead of the benchmark U(2) model we considered its U(iV) general- 
ization, we would see that, as we increase the number of colors N, these 
crossovers become exceedingly sharper and transform into genuine phase 
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transitions in the limit N — > oo. The width of the crossover domain scales 
as 1/N. At finite N there is no discontinuity in physical observables (and 
associated breaking of relevant global symmetries) along the dashed lines in 
Fig. 1. Still both the perturbative spectrum and confining strings are dra- 
matically different in the regimes I, II and III. In particular, if we keep the 
quark mass difference Am = and reduce the FI parameter £ we pass from 
the regime with non-Abelian strings and non-Abelian monopole confinement 
to the regime with Abelian strings and Abelian confinement. The low-energy 
spectrum of the theory in, say, domain I is not mapped onto that in domain 
III. 

The paper is organized as follows. In Sect. 2 we describe our bulk theory 
at large £. In Sect. 3 we make some preliminary remarks on the behavior 
of the theory at small £. In Sect. 4 we briefly review non-Abelian strings 
and discuss the order parameter which can separate a non-Abelian string 
from Abelian one. In Sect. 5 we describe our theory in domain III, while in 
Sect. 6 we consider the limit N — > oo. Section 7 summarizes our findings. 
In Appendix we present in more detail the CP(N — 1) model with twisted 
masses and Z 2 n global symmetry. 

2 The bulk theory: large values of the 
FI parameter 

In this section for convenience of the reader we will briefly outline some basic 
features of the Af = 2 bulk theory we will work with. Since these features 
are general, we will assume the gauge group to be U(iV), and only later we 
will set N = 2. 

We introduce N flavors (each of which is described by two M = 1 super- 
fields, Q and Q). The field content is as follows. The M = 2 vector multiplet 
consists of the U(l) gauge field and the SU(iV) gauge field A a , where 
a = 1, N 2 — 1, and their Weyl fermion superpartners plus complex scalar 
fields a, and a a . The latter are in the adjoint representation of SU(iV). 

The quark multiplets of the SU(iV) xU(l) theory consist of the complex 
scalar fields q and q~Ak (squarks) and their fermion superpartners, all in the 
fundamental representation of the SU(iV) gauge group. Here k = 1,...,N 
is the color index while A is the flavor index, A — 1, ...,N. Note that the 
scalars q kA and q kA form a doublet under the action of the global SU(2) R 
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group. 

The bosonic part of the bulk theory has the form [9] (see also the review 
paper [7]) 



/ 
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+ |V M /| 2 + + V(q A ,q A ,a a ,a 



~A\2 



~A ~ a 



(2.1) 



Here is the covariant derivative in the adjoint representation of SU(iV), 
while 



V» = d,--A,-iA«T°. 



(2.2) 



We suppress the color SU(iV) indices. The normalization of the SU(iV) gen- 
erators T a is as follows 

Tr (T a T b ) = \5 ab . 

The coupling constants g\ and gi correspond to the U(l) and SU(iV) sectors, 
respectively. With our conventions, the U(l) charges of the fundamental 
matter fields are ±1/2. 

The potential V(q A , q A) a a , a) in the action (2.1) is the sum of D and F 
terms, 



V(q A ,q A ,a a ,a) 



2 /1 _ x 2 

2 f abc a b a c + q A T a q A -q A T a ? 



!h_ [ _L fat" t.<> ■ 
2 U 



+ 9 f {ug A - g A r - my 



+ 2g \q A l q + — \ q A q 



1 N fi 

+ - \ K° + v/ 2 m ^ + 2T a a a )g" 

2 a=i ^ 



+ 



> + V2m A + 2T a a a )g 



(2.3) 



Here / afec denote the structure constants of the SU(iV) group, m A is the mass 
term of the A-th flavor, and the sum over the repeated flavor indices A is 
implied. 
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We introduced the FI D-term for the U(l) gauge factor with the FI 
parameter £. 

Now we briefly review the vacuum structure and the excitation spectrum 
in the bulk theory. The vacua of the theory (2.1) are determined by the 
zeros of the potential (2.3). The adjoint fields develop the following vacuum 
expectation values (VEVs): 

<G* + ™)>H5(v:::^)- (24) 

For generic values of the quark masses, the SU(iV) subgroup of the gauge 
group is broken down to U(1) JV_1 . However, in the special limit 

mi = m 2 = ... = ttl 7v , (2.5) 

the SU(iV)xU(l) gauge group remains unbroken by the adjoint field. In this 
limit the theory acquires the global flavor SU(iV) symmetry. 

We can exploit gauge rotations to make all squark VEVs real. Then in 
the case at hand they take the color-flavor locked form 

(q kA ) = ••■ ■■■ ■■■ j , (^) = 0, 

k = 1,...,N A=1,...,N, (2.6) 

where we write down the quark fields as an N x N matrix in the color and 
flavor indices. This particular form of the squark condensates is dictated 
by the third line in Eq. (2.3). Note that the squark fields stabilize at non- 
vanishing values entirely due to the U(l) factor represented by the second 
term in the third line. 

The vacuum field (2.6) results in the spontaneous breaking of both gauge 
and flavor SU(iV) symmetries. A diagonal global SU(iV) c+i ? survives, how- 
ever, 

U(iV) gauge x SU(V) flavor - SV(N) C+F . (2.7) 

Thus, a color-flavor locking takes place in the vacuum. The presence of the 
global S\J(N)c+f group is a key reason for the formation of non-Abelian 
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strings. For generic quark masses the global symmetry (2.7) is broken down 
toU(l)( iV - 1 ). 

Let us move on to the issue of the excitation spectrum in this vacuum 
[16, 9]. The mass matrix for the gauge fields (A^, A^) can be read off from the 
quark kinetic terms in (2.1). It shows that all SU(iV) gauge bosons become 
massive, with one and the same mass 

m w = g 2 \fl. (2.8) 

The equality of the masses is no accident. It is a consequence of the unbroken 
SU(iV)c+F symmetry (2.7). The mass of the U(l) gauge boson is 



rrh = 9i)/-2^ (2-9) 

Thus, the bulk theory is fully Higgsed. The mass spectrum of the adjoint 
scalar excitations is the same as the one for the gauge bosons. This is enforced 
by N = 2 . 

The mass spectrum of the quark excitations can be read off from the 
potential (2.3). We have 4N 2 real degrees of freedom of quark scalars q and 
q. Out of those N 2 are eaten up by the Higgs mechanism. The remaining 
3N 2 states split in three plus 3(iV 2 — 1) states with masses (2.9) and (2.8), 
respectively. Combining these states with the massive gauge bosons and the 
adjoint scalar states we get [16, 9] one long M = 2 BPS multiplet (eight real 
bosonic plus eight fermionic degrees of freedom) with mass (2.9) and iV 2 — 1 
long Af = 2 BPS multiplets with mass (2.8). Note that these supermulti- 
plets come in representations of the unbroken SXJ(N)c+f group, namely, the 
singlet and adjoint representations. 

Now let us have a closer look at quantum effects in the theory (2.1). The 
SU(iV) sector is asymptotically free. The running of the corresponding gauge 
coupling, if not interrupted, would drag the theory into the strong coupling 
regime. This would invalidate our quasiclassical analysis. Moreover, strong 
coupling effects on the Coulomb branch would break SU(iV) gauge subgroup 
down to U(l) iV ~ 1 by virtue of the Seiberg-Witten mechanism [1]. No non- 
Abelian strings would emerge. 

The semiclassical analysis above is valid if the FI parameter £ is large, 

£»A, (2.10) 

where A is the scale of the SU(iV) gauge theory. This condition ensures weak 
coupling in the SU(iV) sector because the SU(iV) gauge coupling does not 
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run below the scale of the quark VEVs which is determined by £. More 
explicitly, 

8n2 7Vln^»l. (2.11) 



gm a 

3 Towards smaller £ 

Below we will see that if we pass to small £ along the line Am = 0, into the 
strong coupling domain, where the condition (2.10) is not met, the theory 
undergoes a crossover transition into the Seiberg-Witten Abelian confine- 
ment regime. In this regime the low-energy perturbative sector contains no 
nontrivial representations of the unbroken SXJ(N)c+f group. Moreover, no 
non-Abelian strings develop. 

The main tool which allows us to identify this crossover transition is 
the presence of the global unbroken SU(iV)(7 + ^ symmetry in the theory at 
hand. First, we note that it is not spontaneously broken in the bulk. If 
it were broken this would imply the presence of massless Goldstone states. 
However, we showed above that the perturbative sector of the theory at large 
£ has a large mass gap of the order of g 2 \Z£, and masses of no states can be 
shifted to zero by small quantum corrections of the order of A. Nor do we 
see the adjoint multiplet of massless Goldstones at small £. 

The presence of the global unbroken SU(N) c +f symmetry means that all 
multiplets should come in representations of this group. We showed above 
that at large £ this is the case indeed: all light states come in adjoint and 
singlet representations ((N 2 — 1) + 1). We will see later that at small £ (in 
the Seiberg-Witten regime) the low-energy spectrum is very different. It 
contains only N states which do not fill any nontrivial representations of 
SXJ(N) C+F . They are all singlets. 

To elucidate the point let us note the following. All (iV 2 — 1) states of the 
adjoint gauge-boson multiplet of SU(iV) C 7 + ^ have degenerate masses (2.8) at 
large £. The presence of the unbroken global S\J(N) c +f ensures that they 
are not split. Imagine that these states were split with small splittings of the 
order of A. Then in the limit of small ^, ^ < A, some of these states could, 
in principle, evolve into (N — 1) light Abelian states while other members of 
the multiplet could acquire large masses, of the order of A (i.e. the photons 
could become light, while the W bosons could become heavy). We stress 
that this does not happen in the theory at hand. The adjoint multiplet is 
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not split at large £ and therefore can disappear from the low-energy spectrum 
at small £ only as a whole. Hence, the light photons in the Seiberg-Witten 
regime at small £ have nothing to do with the diagonal (Cartan) entries of 
the gauge adjoint SXJ(N)c+f multiplet at large £. Similarly, the light dyons 
in the Seiberg-Witten regime at small £ have nothing to do with the light 
quarks q kA , qAk of the non-Abelian confinement regime at large £. The latter 
fill the adjoint representation of SU(iV) c+ ^ (see the discussion above), while 
the former are singlets. 

In order to see what happens with the low-energy spectrum of the bulk 
theory as we reduce £ we use the following method. First we introduce quark 
mass differences (m,A — ms) an d take them large, \m,A — m B \ ^> A. Then we 
can reduce the parameter £ keeping the theory at weak coupling and under 
control. Next, we reach the Coulomb branch at zero £ and use the exact 
solution of the theory [1, 15] to go back to the desired limit of degenerate 
quark masses (2.5). Thus, our routing is: Domain I — > Domain II — > Domain 
III. In domain II the global S\]{N) c +f is lost, and a level crossing occurs. 

The program outlined above will be carried out in full in Sect. 5. 

To conclude this section we briefly review the theory (2.1) at non-zero 
quark mass differences (tua ~ tub) ^ 0, see [10, 7]. At non- vanishing 
(rriA — rn B ) the global S\J (N) c +f is explicitly broken down to U(l)^ jV_1 ^. 
The adjoint multiplet is split. The diagonal entries (photons and their 
Af = 2 quark superpartners) have masses given in (2.8), while the off-diagonal 
states (W bosons and the off-diagonal entries of the quark matrix q kA ) ac- 
quire additional contributions to their masses proportional to (m,A — itib)- 
As we make the mass differences larger, the W bosons become exceedingly 
heavier, decouple from the low-energy spectrum, and we are left with N 
photon states and N diagonal elements of the quark matrix. The low-energy 
spectrum becomes Abelian. 

4 Non-Abelian strings at large £ 

Here we will study the passage from Domain I — > Domain II. At first, we will 
briefly review non-Abelian strings [8, 9, 10, 11] in the theory (2.1), see [7] for 
details. The Abelian Zjv-string solutions break the SU (N)c+f global group. 
Therefore strings have orientational zero modes, associated with rotations 
of their color flux inside the non-Abelian SU(AQ. This makes these strings 
non-Abelian. The global group is broken by the string solution down to 



9 



SU(iV — 1) x U(l). Therefore the moduli space of the non-Abelian string is 
described by the coset 

SV(N) - CP(JV-l). (4.1) 



SU(JV- 1) x U(l) 

The CP(iV — 1) space can be parametrized by a complex vector n l in the 
fundamental representation of SU(iV) subject to the constraint 

n*n l = 1 , (4.2) 

where I = 1, ...,N. As we will show below, one U(l) phase will be gauged 
away in the effective sigma model. This gives the correct number of degrees 
of freedom, namely, 2 (N — 1) . 

With this parametrization the elementary string solution (with the lowest 
winding number) can be written as [10, 5] 

q = -^{(N-l)fa + fa] + (fa-fa) (n -n*-± 



,SU(JV) 



n - n * -\i ) £ ij ^ fNA(r) , 



N J ~ lJ r 2 

A UW = jfS^fir), q kA = 0, (4.3) 



where i — 1,2 labels coordinates in the plane orthogonal to the string axis and 
r and a are the polar coordinates in this plane. For brevity we suppress all 
SU(iV) indices. The profile functions 0i(r) and fair) determine the profiles 
of the scalar fields, while fNA{ r ) and f(r) determine the SU(iV) and U(l) 
gauge fields of the string solution, respectively. These functions satisfy the 
first-order equations [9] which can be solved numerically. 
The tension of the elementary string is given by 

T = 2ttC (4.4) 

Making the moduli vector n l a slowly varying function of the string world 
sheet coordinates [k = 0, 3), we can derive the effective low energy-theory 
on the string world sheet [9, 10, 5]. From the topological reasoning above 
(see (4.1)) it is clear that we will get two-dimensional CP(iV — 1) model. The 
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M = 2 supersymmetric CP(iV — 1) model can be understood as a strong- 
coupling limit of a U(l) gauge theory [17]. Then the bosonic part of the 
action takes the form 

Scp(at-i) = J d 2 x 1 2(5 \V k n l \ 2 + -^F 2 kl + ^\d k a\ 2 

+ Af3\a\ 2 \n l \ 2 + 2e 2 !3 2 {\n l \ 2 -l) 2 }, (4.5) 

where V k = d k — %A k while a is a complex scalar field. The condition (4.2) 
is implemented in the limit e 2 — > oo. Moreover, in this limit the gauge 
field A k and its M = 2 bosonic superpartner a become auxiliary and can be 
eliminated by virtue of the equations of motion, 

A k = ~n* e Sk n e , a = 0. (4.6) 

The two-dimensional coupling constant (3 here is determined by the four- 
dimensional non-Abelian coupling via the relation 

/> = £■ («) 

The above relation between the four-dimensional and two-dimensional cou- 
pling constants (4.7) is obtained at the classical level [9, 10]. In quantum 
theory both couplings run. In particular, the CP (N — 1) model is asymp- 
totically free [18] and develops its own scale A CT . The ultraviolet cut-off of 
the sigma model on the string worldsheet is determined by g2\^- Equa- 
tion (4.7) relating the two- and four-dimensional couplings is valid at this 
scale, implying 

A? = g?t*e- S £=A' f . (4.8) 

Note that in the bulk theory per se, because of the VEVs of the squark 
fields, the coupling constant is frozen at (y^v^j there are no logarithms below 
this scale. The logarithms of the string worldsheet theory take over. More- 
over, the dynamical scales of the bulk and worldsheet theories turn out to be 
the same [10]. 

The CP(N - 1) model was solved by Witten in the large-A limit [19]. 
We will briefly summarize Witten's results and translate them in terms of 
strings in four dimensions [10]. 
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Classically the field n can have arbitrary direction; therefore, one might 
naively expect a spontaneous breaking of SU(iV) and the occurrence of mass- 
less Goldstone modes on the string world sheet. Well, this cannot happen 
in two dimensions. Quantum effects restore the symmetry. Moreover, the 
condition (4.2) gets in effect relaxed. Due to strong coupling we have more 
degrees of freedom than in the original Lagrangian, namely all N fields n 
become dynamical and acquire masses A a . 

Deep in the quantum non-Abelian regime the CP(iV — l)-model strings 
carry no average SU(iV) magnetic flux. To see that this is indeed the case, 
note that the SU(iV) magnetic flux of the non-Abelian string (4.3) is given 
by 

J d 2 x (F 3 *) SU(JV) = 2tt (n ■ n* - 1) , (4.9) 

where 

F; = \e ijk F jk , (i,j,k = 1,2,3). (4.10) 

As was shown by Witten [19], in the CP(iV — 1) model strong quantum 
fluctuations of n e result in 

(n £ )=0, (4.11) 

implying, in turn, that the average SU(iV) magnetic flux of the non-Abelian 
string vanishes. We will use this circumstance later, to distinguish between 
large-£ non-Abelian and small-£ Abelian Z N strings in the Seiberg-Witten 
regime below the crossover. The latter do carry the magnetic flux directed 
inside the Cartan subalgebra of SU(iV). The CP(N — 1) model has iV vacua 
[19]. They are interpreted in the problem at hand as N different elementary 
non-Abelian strings. These iV vacua differ from each other by the expectation 
value of the chiral bifermion operator, see e.g. [20]. At strong coupling the 
chiral condensate is the order parameter for Z^ breaking (instead of the 
flux, see Appendix). The U(l) chiral symmetry of the CP(iV — 1) model is 
explicitly broken to a discrete Z 2 n symmetry by the chiral anomaly (for a 
discussion of the global symmetry on the string world sheet see Appendix). 
The bifermion condensate breaks Z 2 n down to Z 2 . That's the origin of the 
A"-fold degeneracy of the vacuum state. 

Now, to make our consideration simpler we will focus on the simplest 
case N = 2. For arbitrary N the emerging dynamical pattern is similar. The 



12 



solution for the non-Abelian string (4.3) in the N = 2 case takes the form 
q = ^(0i + 2 ) + ^ a (0i-0 2 ), 9 = 0, 

(4.12) 

where S a (a = 1, 2, 3) is a real moduli vector subject to the constraint 



At(x) = S a e tJ ^ f NA {r) , A(x) = e i3 ^ /(r) , 



{S a f = 1 . 

Its relation to the complex vector n e is as follows 

S a = fi r a n. 



(4.13) 



(4.14) 



We have CP(1) as the effective world-sheet theory. It is equivalent to the 
0(3) sigma model. In terms of real vector S a the bosonic part of world-sheet 
theory has the form 



s = p 



J d 2 x±(d k s a y 



(4.15) 



Now let us introduce quark mass differences (itla -ffls). In the N = 2 
case we have just one (generally speaking complex) parameter 



Am = mi — m 2 . 

The vacuum expectation values of the adjoint field reduce to 

Am , , nz mi + m 2 



V2' 



(a) = -y/2 



(4.16) 



(4.17) 



see (2.4). The non-Abelian string (4.12) is no longer a solution of the first- 
order equations for arbitrary S a . The global SU (2) c+f is explicitly broken 
down to U(l) by Am ^ 0. Nevertheless, if we keep Am small, we can 
consider S a as quasimoduli, with a shallow potential on the CP(1) moduli 
space. The string solution (4.12) in this case should be supplemented by a 
nontrivial profile for the adjoint field [10, 7], 



a a = 



Am 



+ 5"* 5-3(1- fl ) 

<P2 <P2 



01 



(4.18) 
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Plugging this modified string solution in the action of the theory gives [10, 7] 
the effective string world-sheet theory: M = 2 CP(1) model with twisted 
mass [21]. The bosonic part of the action is 

Scp(i) =fij d*x {1 (d k S°) 2 + ^ (1 - Si) } . (4.19) 

This is the only functional form that allows Af = 2 completion. The mass- 
splitting parameter Am of the bulk theory exactly coincides with the twisted 
mass of the world-sheet model. 

The CP(1) model (4.19) has two vacua located at S a = (0, 0, ±1). Clearly 
these two vacua correspond to two elementary Z 2 strings. 

With non-vanishing Am we can introduce a gauge invariant quantity 
which measures the SU(2) non-Abelian flux of the string. We define 2 

$ = J d 2 xa a F; a . (4.20) 

This order parameter will be used below to distinguish between different 
regimes of the theory. 

Substituting (4.18) and (4.9) into (4.20) we get 

$ = _ 27r ^ 1 S 3 (4.21) 
v2 

At small Am, Am <C A CT , the fields S a strongly fluctuate and (S* 3 ) = (see 
(4.11)). Therefore, 

-> at Am < A,j , (4.22) 

where the subscript I refers to the non-Abelian domain at large £ and small 
Am, as it is indicated in Fig. 1. 

Instead, at large Am (Am 3> A CT ) the 0(3) sigma model (4.19) is at weak 
coupling. Fluctuations are small, and the S a acquires vacuum values at the 
north and south poles of the £2 sphere, (S a ) = (0, 0, ±1). As a result 

A 777 

(<&)n-=F27r-^, (4.23) 
where the subscript II marks domain II in Fig. 1. 



2 Thc subscript 3 indicates the direction along the string axis, cf. Eq. (4.9). 
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Figure 2: Flux (4.20) as a function of Am in domains I and II. 



We see that the behavior of the string flux $ drastically changes as we 
pass from the non-Abelian domain I of large £ and small Am to the Abelian 
domain II of large Am, see Fig. 2. At large Am this theory is in the weak 
coupling regime and fluctuations are small, (S 3 ) ±1 and the flux $ is 
given by (4.23). At small Am the world sheet theory is in the strong coupled 
quantum regime, fluctuations are large and the vector S a is smeared over 
the whole sphere. Therefore, (S 3 ) ~ and (<&) ~ 0. The crossover between 
these two regimes is at Am ~ A a . Note, that the drastic change of behavior 
in the world-sheet CP(1) model is correlated with the dynamics of the bulk 
theory. 

In Sect. 2 we saw that the perturbative spectrum of the bulk theory is 
different in these two domains: it is essentially non-Abelian in the domain 
I while in the domain II the W bosons become exceedingly heavier and 
decouple from the low-energy spectrum: we are left with N photons and 
their superpartners, the diagonal elements of the quark matrix. The pattern 
repeats itself at the nonperturbative level: the non-Abelian strings evolve 
into the Abelian strings as we increase |Am|. 

Later we will see that the crossover becomes exceedingly more pronounced 
as we increase the number of colors N. In the limit N — > oo the crossover 
evolves into a genuine phase transition. Note also that in nonsupersymmetric 
theories we do have a phase transition between the phase with non-Abelian 
strings at small \m,A — itib\ and the phase with the Abelian strings at large 
\rriA — m B \ [6]. It is related to the restoration of the broken discrete Z N 
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symmetry at small \m,A — tub]- I n the supersymmetric theory at hand the 
discrete symmetry is always broken (by VEVs of n at large — m#| 
or by the two-dimensional bifermion condensates at small \rriA — ttlb\, see 
Appendix). Therefore in the supersymmetric case we have a crossover rather 
than a phase transition. 

To conclude this section, we briefly review the world-sheet theory on 
the non-Abelian string for generic N. It is described by the twisted-mass- 
deformed CP(N — 1) model. It can be nicely written [22] as a strong coupling 
limit of a U(l) gauge theory. With twisted masses of the n e fields taken into 
account, the bosonic part of the action (4.5) takes the form 



2(3\V k n e \ 2 + 



+ A{3 



rru 
V2 



\n 



£\2 



+ 2e 2 2 (\n e \ 2 -l 



(4.24) 



where 



rri£ = mi — m, m 



N ^ 



(4.25) 



and the sum over £ in (4.24) is implied. 



5 The theory in domain III 

Now, we will consider the passage from domain II to domain III. In order 
to study the theory in the regime III (see Fig. 1) we first assume the quark 
mass differences to be large. Then the theory stays at weak coupling and we 
can safely decrease the value of the FI parameter £. Next, we use the exact 
Seiberg-Witten solution of the theory on the Coulomb branch (at £ = 0) to 
pass from regime II to regime III. To simplify our discussion, we will consider 
here only the case N = 2. 

3 At N > 2 the above discrete symmetry of the Lagrangian takes place under a special 
choice of the mass parameters, see Appendix. 
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5.1 The r = 2 quark vacuum 

Our first task is to identify the r = 2 quark vacuum (which we described 
semiclassically above) using the exact Seiberg-Witten solution [1, 15] 4 . The 
Seiberg-Witten curve for the U(2) gauge theory with Nf — 2 flavors has the 
form 

y 2 = (x- fc) 2 (* - 2 ) 2 - 4A 2 (x + ^ )(x + (5.1) 



where 0i and 2 are gauge-invariant parameters on the Coulomb branch. 
Semiclassically 

0i ~ ai = i(a + a 3 ), 2 ~ a 2 = ^(a - a 3 ). (5.2) 
Let us make a shift in the variable x introducing a new variable z, 

Tfl 1 

x = i= + z, m— -(m 1 + m 2 ). (5.3) 

V 2 2 

With Am > A we identify the r = 2 singularity, the point where both 
quarks g 11 and g 22 become massless. Upon switching on £ 7^ 0, this r = 2 
singularity turns into the r = 2 vacuum we considered in the semiclassical 
approximation in the previous sections. 

It turns out that in the r = 2 vacuum 

Tfl 

0! + 2 = -2-j= . (5.4) 

We parametrize the deviations of X and 2 from their mean value — ^= by 
a new parameter 0, 

m 

0! = "^ + 0, 

02 = ( 5 - 5 ) 



4 Thcsc solutions were obtained by Seiberg and Witten in the SU(2) gauge theories. 
Generalizations to SU(A r ) were obtained in [23, 24, 25, 26]. 
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With this parametrization the curve (5.1) reduces to ! 



(z - <p)\z + 0) 2 - 4A 2 (z + 



Am 
2^2 



2\2 



= (^-0 2 ) 



4A 2 



Am 2 



(5.6) 



Next, we look for the values of the parameter which ensure that this 
curve has two double roots associated with two quarks being massless. This 
curve is a perfect square 

2 



y 2 = 



1 /Am 2 
4 V 2 



+ 4A' 



see [22], at 



, 1 /Am 2 
0=-^— -4A 2 



(5.7) 



(5.8) 



In fact, there are two solutions with plus and minus signs in front of the 
square root above. They correspond to 0i and 02, namely 



ml/ Am 2 

0i = f= \ 4A 2 , 

r V2 2 V 2 



02 = 



m 1 / Am 2 

■— f= + -\ 4A 2 . 

y/2 2V 2 



(5.9) 



In the semiclassical limit Am ^> A these solutions reduce to 

X m l A 1712 / r ia\ 

0x^-—, 2 (5.10) 

which coincides with Eq. (2.4). This means that we correctly identified the 
r = 2 quark vacuum where two quarks q 11 and g 22 are massless, see Sect. 2. 
Two double roots of the curve in the quark vacuum are 



ei = e 2 = 



e 3 = e 4 = 



m 
71 



Am 2 



+ 4A 2 , 



m 1 Am 2 

+ -\ h4A 2 . 

y/2 2 V 2 



(5.11) 



5 Am 2 is a shorthand for (Am) 2 . 



18 



The Seiberg-Witten exact solution of the theory relates VEVs of the fields 
a, a 3 and ap, a 3 D (which, in turn, determine the spectrum of the BPS states 
on the Coulomb branch) to certain integrals along a and (3 contours in the 
x-plane [1, 15]. Say, the derivatives of VEVs of a\ and a 2 are given by the 
following integrals along the a contours: 



da i 


1 


90i 


2tH 


da i 


1 


902 




da 2 


1 


90i 




da 2 


1 


902 


2tH 



IE -02 

ax 



y 

X-(f> 

ax 

y 



x-<j> 
ax — 



2 



y 

dx^±, (5.12) 

v 

while the derivatives of ao's are given by similar integrals along the (3 con- 
tours. 

The presence of two massless quarks q 11 and q 22 in the r = 2 vacuum at 
Am ^> A implies 

ai + ^ = 0, «2 + ^= = 0. (5.13) 

Thus, the fields ai, 02 are regular at the singularity while the fields an 
have logarithmic divergences related to the /3 functions of the low-energy 
U(l)xU(l) theory. This ensures that the a\ contour should go around the 
roots ei, e 2 while the a 2 contour should go around the roots e^, e±. In the 
r = 2 vacuum (5.11) both contours shrink and produce regular a's. The basis 
of the a and (3 contours is shown in Fig. 3. Here we consider our U(2) theory 
as a two-flavor SU(3) gauge theory broken down to U(2) at a very high scale. 
In terms of the SEiberg-Witten curve this corresponds to extra two roots of 
the SU(3) curve being far away from four roots of the U(2) curve (5.1). 

As a double check of our identification of the quark vacuum let us cal- 
culate the derivatives (5.12) in the semiclassical approximation Am ^> A. 
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Figure 3: Basis of the a and contours of our U(2) gauge theory viewed as an 
SU(3) theory broken down to U(2). Two extra roots of the SU(3) theory are far 
away in the x plane. 



Substituting (5.7) into (5.12) we get at Am 3> A 



da i 

Wi 

da 2 
Wi 







da i 
d<p 2 

da 2 

d(j) 2 



0, 



1. 



(5.14) 



This is in accord with (5.2) and confirms our choice of the a contours in 
Fig. 3. 

To conclude this subsection we note that the monopole singularity (the 
point on the Coulomb branch where the SU(2) monopole becomes massless, 
a 3 D = 0) corresponds to shrinking of the (fli—fa) contour, i.e, in other words, 
to ei = e 3 . 



5.2 Monodromies 

Let us study how the quantum numbers of massless quarks q 11 and q 22 change 
as we reduce | Am| and go from domain II into domain III where the theory is 
at strong coupling. The quantum numbers change due to monodromies with 
respect to Am. The complex plane of Am has cuts and when we cross these 
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cuts, the a and ap fields acquire monodromies and the quantum numbers of 
states change accordingly. Monodromies with respect to quark masses were 
studied in [27] in the theory with the SU(2) gauge group using a monodromy 
matrix approach. 

Here we will investigate the monodromies in the U(2) theory with two 
quark flavors using a slightly different approach, similar to that of Ref. [28]. 
If two roots of the Seiberg-Witten curve coincide, the contour which goes 
around these roots shrinks and produces a regular potential. Say, as was 
discussed above, at Am > A we have two double roots e\ = e 2 and = e± 
in the r = 2 vacuum. Thus, two contours a\ and a 2 shrink (see Fig. 3), and 
potentials a\ and a 2 are regular. This is associated with masslessness of two 
quarks, see (5.13). 

Instead, in the monopole singularity e\ = e^, thus the (f3\ — f3 2 ) contour 
shrinks producing a regular a 3 D . This is associated with the masslessness of 
the SU(2) monopole, a 3 D = [1]. 

If we decrease \Am\ and cross the cuts in the Am plane, the root pairing in 
the given vacuum may change. This would mean that a different combination 
of a and becomes regular implying a change of the quantum numbers of 
the massless states in the given vacuum. To see how it works for our r = 2 
vacuum we go to the Argyres-Douglas (AD) point point [29, 30]. The AD 
point is a particular value of the quark mass parameters where more mutually 
nonlocal states become massless. In fact, we will study the collision of the 
r = 2 quark vacuum with the monopole singularity. We approach the AD 
point from domain II at large Am. We will show below that as we pass 
through the AD point the root pairings change in the r = 2 vacuum implying 
a change of the quantum numbers of the massless states. Two massless quarks 
transform into two massless dyons. 

To be more precise, we collide the r = 2 vacuum with two massless quarks 
with the quantum numbers 

(n e ,n m ;n 3 e ,n 3 J = (1/2, 0; 1/2, 0), 

(n e ,n m ,n 3 e ,n 3 J = (1/2, 0; -1/2, 0) (5.15) 

with the monopole singularity with 

(n e ,n m ,n 3 e ,n 3 m ) = (0,0; 0,1) (5.16) 

where the monopole becomes massless. Here n e and n m denote electric and 
magnetic charges of a state with respect to U(l) gauge group, while n 3 and r? m 
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stands for electric and magnetic charges with respect to the Cartan generator 
of SU(2) gauge group (broken down to U(l) by Am). 

As was already mentioned, the (0, 0; 0, 1) monopole is massless if e\ = e^. 
Equation (5.11) shows that this can happen in the r = 2 vacuum only if all 
four roots of the U(2) curve coincide at 

Am 2 = -8A 2 , ei = e 2 = e 3 = e 4 = -= ■ (5.17) 

v2 

This is the position of the AD point where both, the quarks and the SU(2) 
monopole become simultaneously massless. 

In order to see how the root pairings in the r = 2 vacuum change as 
we decrease |Am| and pass from domain II into domain III through the AD 
point (5.17) we have to slightly split the roots by shifting <fi from its r = 2 
solution (5.8). Let us take 

^i(^ 2 ) + >' (5 ' 18) 

where 5 is a small deviation. Then the curve (5.1) can be approximately (at 
small 5 ) written as 



y ~ 



4 

Now all four roots split as follows: 

ei 



2 1 /Am 2 ax2 
z — - I — — + 4A 2 



2 

2 



5 2 . (5.19) 



m i—z r m j—z 

= -= + V/i 2 + $ , e 2 = -= + V/i 2 - 



777 , 777 , 

e 3 = + e 4 = --^-v / ^ r 5 ) (5-20) 

where we introduced a shorthand notation 



"4v^ +4A2 - (5 - 2i) 

This parameter vanishes at the AD point. 

In order to pass through the AD point from domain II into domain III 
we decrease |Am| keeping Am pure imaginary, 

Am=\Am\e i *. (5.22) 
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Figure 4: As we decrease \Am\ (keeping Am imaginary) and pass through the 
AD point, the roots 61,2,3,4 move in the x plane. 

Then \x goes along the imaginary axis towards the origin (which is the AD 
point) and below the AD point increases along the positive axis. We also fix 
the parameter 8 to be imaginary too, 5 = \S\ e*2 . This is convenient as all 
four roots stay split at any \Am\. 

As we decrease \Am\ the roots (5.20) move as shown in Fig. 4. We see 
that the root pairings in the r = 2 vacuum change. Namely, at large |Am| 
we have (at 5 = 0) 

ei = e 2 , e 3 = e 4 , (5.23) 

which, as was explained above, corresponds to shrinking of the a± and a 2 
contours and masslessness of two quarks (5.15). Below the AD point at small 
I Am I we have 

e 2 = e 3 , ei = e 4 , (5.24) 
which corresponds to shrinking of the contours 

A - fa + 011 -> 0, -ft + fa + a 2 -> . (5.25) 

This means that massless quarks in the r = 2 vacuum transformed into 
massless dyons D\ and D 2 with the quantum numbers 

D l : (1/2,0; 1/2,1), D 2 : (1/2, 0; -1/2, -1) . (5.26) 
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We see that the quantum numbers of the massless quarks in the r = 2 
vacuum after the collision with the monopole singularity get shifted, the shift 
being equal to ±(monopole magnetic charge). 

The monodromy discussed above implies 

ai -> ai + a 3 D , a 2 -> a 2 - a 3 D , a 3 -> a 3 + 2a|, . (5.27) 

Therefore, the conditions (5.13) for masslessness of the q 11 and q 22 quarks 
are replaced in domain III by the conditions of masslessness of the dyons Di 
and D 2 , namely, 

, m-i , mo 

a 1 + a 3 D + —L = 0, a 2 -a 3 D + -^ = 0. (5.28) 
v2 v 2 

5.3 The low-energy theory 

In this subsection we present the low-energy theory in the r = 2 vacuum in 
domain III at small £ and small |Am| (below the AD point). It should be 
stressed that none of the fields in this low-energy theory belong to nontrivial 
representations of S\J(2) c +f- 

As we already know, the massless quarks g n and q 22 transform into the 
massless dyons D\ and D 2 . The latter interact with two photons. According 
to the dyon quantum numbers (5.26) one of these photons is 

A„ (5.29) 

while the other photon is the following linear combination: 

B » = -j$ (4 + 2A T) ■ ( 5 - 3 °) 

In fact, these are the only light states to be included in the low-energy 
effective theory in domain III. All other states are either heavy (with masses 
of the order of A) or decay on curves of marginal stability. In the case at 
hand CMS is located around the origin in the Am complex plane and goes 
through the AD point [31]. In fact, the W bosons of the underlying non- 
Abelian gauge theory, as well as the off-diagonal states of the quark matrix 
q kA , decay on CMS. Let us illustrate this statement, say, for the W bosons. 
To this end we can go to the AD point. At this point we have for the W-boson 
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mass 



m w = V2\a 3 \ = V2 I (a 3 + a 3 D ) - a 3 D \ 

= V2(\a 3 + a 3 D \ + \a 3 D \) = m M + m D , 



(5.31) 



where uim and mo are the masses of the SU(2) monopole and SU(2) dyon 
with charges (0,0; 0,1) and (0,0; 1,1), respectively. This relation is valid 
at the AD point just because the monopole becomes massless at this point, 
a 3 D = 0. It means that the W-boson decays into the SU(2) monopole and 
dyon at this point and is not present in domain III, in full accordance with 
the analysis of the SU(2) theory in [27]. 

Taking this into account we write the effective low-energy action of the 
theory in domain III as follows: 



■>iii 



+ 



+ 



K D i\ 



v 2 A 



V(D,D,b,a)] , 



where 



b = -j={a 3 + 2a 3 D ) 



(5.32) 



(5.33) 



is the scalar Af = 2 superpartner of the photon (5.30) while F^ u is the field 
strength of the U(l) gauge field B^. Covariant derivatives are defined in 
accordance with the charges of the dyons D 1 and D 2 . Namely, 



K = d,-r{^A,-\Al-Af\=d,-\{A,-^B,).^.U) 



The coupling constants g 1 and g 2 correspond to two U(l) gauge groups. The 
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potential V(D, D, b, a) in the action (5.32) is 

5~g 2 < 



V(D,D,b,a) 



DA - LDi - \D 2 \ + \D 2 



+ ?r (l^il 2 - \Di\ 2 + \D 2 \ 2 - \D 2 12 



+ 



5~9 2 2 



D 1 D 1 - D 2 D 2 



2 gl 



D X D X + D 2 D 2 



+ X - | a + VEb + yfirm 2 (|£>i| 2 + I^Dil 2 ) 

+ a- v^ + v^ 2 (| J D 2 | 2 + l^l^j . (5.35) 

Now we are ready move to the desired limit of equal quark masses, Am = 
0. In this limit the global SU(2) C+F symmetry is restored in the underlying 
theory. The vacuum of the theory (5.32) is located at the following values of 
scalars a and b: ^ 

a = —V2m, \fbb = 



V2' 

while the VEVs of dyons are determined by the FI parameter £, 
D^y/Z, D 2 = ^l, D 1 = D 2 = 0. 



(5.36) 



(5.37) 



Thus, the U(l)xU(l) gauge group is broken by dyon condensation. Both, 
photons and dyons become massive, with masses proportional to y/£. In 
particular, at Am = the vacuum value of b vanishes. 

Note also that the theory (5.32) is the Abelian U(l)xU(l) gauge theory 
and hence is not asymptotically free. It stays at weak coupling at small £. 

The low-energy theory (5.32) does not seem to have any global SU(2) 
symmetry. However, the underlying theory does have a global SU(2) sym- 
metry in the limit Am = 0. As was explained in Sect. 2, this global SU(2) 
is not broken in domain I at large £. This symmetry is realized in a color- 
flavor- locked form in this domain (see Eq. (2.7)), and no Goldstone bosons 
are present. We showed that no massless states are present in domain III at 
non-zero £, (and no light states other than two dyons and two photons dis- 
cussed above); therefore, the global SU(2) cannot be spontaneously broken 
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in this domain either. The only way out of this puzzle is to conclude that 
the SU(2) global symmetry is realized trivially in the low-energy description 
(5.32), i.e. that all states in (5.32) are singlets of the unbroken flavor SU(2). 

This means, as was already mentioned in Sect. 3, that the photon 
which appears in domain III has nothing to do with the third component of 
the SU(2) gauge field A® of domain I. At Am = the former is a singlet of 
the global SU(2), while the latter is a component of a triplet. 

Moreover, dyons D\ and D2 present in domain III have nothing to do 
with diagonal entries of the quark matrix q kA of domain I. Dyons are singlets 
while the quarks q kA form the singlet and triplet states. 

Since we have a crossover between domains I and III rather than a phase 
transition, this means that in the full theory triplets become heavy and de- 
couple as we pass from domain I into domain III along the line Am = 0. 
Moreover, some composite singlets, which are heavy and invisible in domain 
I become light in domain III and form dyons Di >2 and photon (level cross- 
ing). Although this crossover is smooth in the full theory, from the stand- 
point of the low-energy description the passage from domain I into domain 
III means a dramatic change: the low-energy theories in these domains are 
completely different, in particular, the degrees of freedom in these theories 
are different (non-Abelian in domain I vs. Abelian in domain III). 

5.4 Strings in domain III 

It is obvious that the low-energy theory (5.32) have Z 2 string solutions in 
the vacuum (5.36), (5.37). Say, the D\ dyon can have a winding at infinity. 
In this case the string solution has the following behavior at r — > 00: 

Di ~ e ia ^, D 2 ~^l, 

Ai ~ di(x, VbBi^dia, (5.38) 

where the indices i — 1,2 denote the plane orthogonal to the string axis and 
r and a are polar coordinates in this plane. Another elementary string can 
be obtained from the one in (5.38) by the replacement Di — > D 2 , D 2 — > D x 
and Bi — > —B{. 

These Z 2 elementary strings are BPS-saturated. Their tensions are given 
by the formula (4.4) in the same way as the tensions of elementary strings in 
domains I and II. 
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The Z 2 strings are Abelian (of the Abrikosov-Nielsen-Olesen type [32]) 
in domain III. They do not have any orientational moduli, in contrast with 
non-Abelian strings in domain I. 

Let us calculate the gauge invariant non-Abelian flux (4.20) for these 
strings. In Abelian domain II at large \Am\ 

a a F; a ^a 3 Ff. (5.39) 

With | Am | decreasing, as we pass through monodromies, we get 



- ~ 3 1 2a 3 D = V5b, (5.40) 



a — > a 



AJ - VhB^. (5.41) 
Equation (5.41) follows from Eq. (5.40) by M = 2 supersymmetry. Therefore 

$ m = J d 2 x(V5b) (V5F; B ) . (5.42) 

Equation (5.36) gives 

Vhb= -Am/V2 

in the r = 2 vacuum, while the flux of the field of the Z 2 strings can be 
read off from Eq. (5.38). In this way we arrive at 

<*>m - T27T ^ . (5.43) 

We see that the string flux in domain III is given by the same formula as 
in domain II. This is a flux of the Abelian string. The non-Abelian part of the 
flux is directed in the Cartan subalgebra of the gauge group. No orientational 
moduli appear. In contrast, in domain I the flux of the non-Abelian string is 
proportional to the orientational vector S a . At small \Am\ the expectation 
value (S a ) — > 0, and the string flux is averaged to zero, see (4.22). Domains 
I and III are separated by a crossover at £ ~ A 2 . 

Let us also mention one more dramatic distinction of nonperturbative 
spectra in domains I and III at Am = 0. The confined SU(2) monopoles 
(with quantum numbers (5.16)) are the junctions of two different elementary 
strings in both domains. In the non-Abelian domain I the confined monopoles 
are seen as kinks of the world sheet CP(iV — 1) model [10, 33, 11]. As was 
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shown by Witten [19], deep in the quantum regime at (m^-mg) = the kink 
of the CP(iV — 1) model is described by the field n l and therefore acquires a 
global flavor quantum number with respect to the unbroken SU (N)c+f- In 
fact, the kink/monopole is in the fundamental representation of this group (a 
doublet in the case N = 2) [19, 34]. Therefore, a meson formed by a monopole 
connected to an antimonopole by two strings (see the review paper [7] for 
details) belongs to the singlet or adjoint representations of the global SU(iV) 
(singlet or triplet of SU(2) for N — 2). 

Clearly, in domain III the monopole confined by strings does not acquire 
global quantum numbers. It is in the singlet representation of SU(2). Hence, 
a meson formed by a monopole connected to an antimonopole by two strings 
is a singlet too. Thus, in the nonperturbative spectra of the theory we 
observe the same phenomenon which was seen in the perturbative spectra: 
triplets of global SU(2) present at low energies in domain I are lifted and do 
not appear in the low-energy description in domain III. Both perturbative 
and nonperturbative states in domain III are singlets of the unbroken global 
SU(2). 

6 The phase transition at AT — > oo 

In this section we will consider the N dependence of the crossover transitions 
(see Fig. 1) in parameters £ and (mA—rn B ). We will show that in the large- N 
limit the crossovers become exceedingly sharper and at N = oo transform 
into genuine phase transition. We will start from the crossover in (m,A — tub) 
at large £ (i.e. the passage from domain I to domain II). 

This crossover in the nonperturbative sector of the theory can be seen 
as a crossover in the effective CP(iV — 1) model (4.24) on the world sheet 
of the non-Abelian string, see Sect. 4. As was already explained, at large 
quark mass differences the CP(iV — 1) model is at weak coupling. The VEV 
of the vector n l does not vanish. If we make a special choice for the mass 
parameters 

m k = m e^~\ k — l,...,N, (6.1) 

where m is a single common parameter (which we will take to be real) our 
theory has a discrete Z 2 n symmetry, see Appendix for further details. In 
fact, (n l ) is an order parameter for the spontaneous breaking of this Z 2N 
symmetry down to Z 2 . 
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At weak coupling, at large m , dynamics can be described as follows. The 
action (4.24) contains a term 



mi 

a — 

V2 



n l \ 2 . (6.2) 



At weak coupling the field n can develop a VEV if a reduces to a particular 
mass parameter, 

*=^=, n l = ^2p S l \ (6.3) 

where k = 1,...,N labels N different vacua (i.e. the elementary Zn strings 
of the bulk theory) and we rescaled the field n l in (4.24) to make its kinetic 
term canonic, namely, n l — > n l /y/2f3. 

As we reduce the value of m the vacuum expectation value of the n l 
field becomes smaller and tends to zero at the left boundary of domain II. 
Simultaneously, the VEV of a is no longer given by the mass, as in Eq. (6.3). 
In fact, o determines the bifermion condensate; \a\ becomes of the order of 
A CT at itlq — > 0. In both limits the symmetry is broken. This is the reason 
why two domains, I and II, are separated by a crossover rather than a phase 
transition. 6 

In order to study the crossover at any N (rather than at N = oo) we can 
use the description of the supersymmetric CP(iV — 1) model in terms of an 
exact superpotential [17, 22]. Upon integrating out n l fields the model can be 
described by an exact twisted superpotential of the Veneziano-Yankielowicz 
type [35] 

where S is a twisted superfield [17] (with a being its lowest scalar compo- 
nent) and we ignore here the 6 dependence (9 stands for the vacuum angle). 
Minimizing this superpotential with respect to a we find 

N 

\J{V2a-m l )=K N a1 (6.5) 



6 In the nonsupersymmctric case the VEV of a vanishes in the domain I, and the 
symmetry is restored [5, 6]. In this case we do have a phase transition between domains I 
and II. 
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where A a is the scale parameter of the CP(iV — 1) sigma model under con- 
sideration. 

Let us examine this equation, determining the VEV of the field a at finite 
rather than infinite N. If iV is fixed, it is readily seen that at large \mi\ (i.e. 
m A CT ) the solution for a coincides with one of the masses, in accordance 
with our semiclassical analysis, see Eq. (6.3). In the opposite limit of zero 
masses (m = 0) 

1 , 2-7T i k . 

a=—A a e—, (6.6) 

where k = 1,...,N marks N distinct vacua. As was mentioned above, the 
Z N symmetry is spontaneously broken at any m . 

As we increase the value of mo, the vacuum expectation of a smoothly 
interpolates between the regime (6.6), where the order parameter which dis- 
tinguishes different vacua of the CP (AT — 1) model (i.e. different elementary 
strings of the bulk theory) is a bifermion condensate ~ a, and the regime 
(6.3) where a is determined by one of the masses mi, while n l develops a 
VEV. For finite N the solution for a is a smooth function of m . Thus, this 
is a crossover that takes place between domains I and II. 

If we increase A" this crossover becomes more pronounced. Let us study 
Eq. (6.5) at large N. To simplify our analysis let us consider A" = 2 P , where 
p is an integer. Then Eq. (6.5) can be rewritten as 

(V2a) N -m% = A?. (6.7) 

This equation has the following perfectly smooth solution: 

a = -Le^* (m» + A^ N . (6.8) 

However, at N — > oo the above function takes the form 



V2 



e~ l x { . (6.9) 

A a , m < A a 



Corrections to this expression are exponential in (-N). 

We see that the solution for a develops a discontinuity in the first deriva- 
tive with respect to m . The crossover becomes a phase transition in the 
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limit N = oo. We stress that this phase transition is an artifact of the large- 
iV approximation and is not related to a change in the pattern of realization 
of any symmetry. 

The solution (6.9) for a ensures the following behavior of the vector n l in 
the iV — > oo limit: 

(n l ) = I (6.10) 
[ 0, m < A a , 

where k = 1, N. The renormalized coupling /3 ren tends to zero at m = A a 
[6]; thus, the VEV of develops a discontinuity in the first derivative with 
respect to m . 

This solution implies that the gauge invariant non-Abelian flux of the 
non-Abelian string strictly vanishes in domain I, 

($)i = (6.11) 

at A^ = oo while in domains II and III it is given by an U(iV) generalization 
of Eq. (5.43). Namely, 

(Q)n = (Q)m = -2nV2m k (6.12) 

for the k-th elementary Z N string, k — 1, ...,iV, see Eqs. (4.9) and (2.4). 

The result (6.11) is exact at A" = oo. Thus, at A^ = oo the string 
flux (4.20) develops a discontinuity as we pass from domain I to domains II 
or III. This implies that both crossovers in £ and m transform into phase 
transitions. 



7 Discussion 

In this paper we considered the r = N vacuum in Af = 2 supersymmetric 
QCD with the U(N) gauge group and Nf = N flavors. We demonstrated that 
this theory exhibits a crossover transition in £, see Fig 1. Namely, at large £ 
in domain I the theory is in the non-Abelian confinement regime, it has A^ 2 
degrees of freedom (gauge bosons and quarks) at low energies and supports 
non-Abelian strings. In contrast, at small £ the theory passes into the Abelian 
Seiberg-Witten regime III. The low-energy effective description includes A^ 
degrees of freedom (dyons and dual photons) and supports Abelian strings. 
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We have shown that non-Abelian gauge bosons and quarks in domain I have 
nothing to do with Abelian dyons and photons of domain III. These states 
belong to different representation of the unbroken global flavor group SU(N). 

Although in this paper we considered a particular vacuum in a specially 
chosen version of M = 2 SQCD (where the global SU(iV) symmetry remains 
unbroken due to the color-flavor locking) we believe that our results are quite 
general. It seems plausible that many Abelian vacua of the Seiberg-Witten 
type in Af = 2 SQCD exhibit crossover transitions into non-Abelian regimes 
as we increase the FI parameter £. Usually we just do not have appropriate 
extra parameters (such as the quark mass differences in our example) which 
would allow us to study these crossovers. 

The lesson is that, generally speaking, non-Abelian strings can smoothly 
evolve into Abelian strings and vice versa. At the same time the correspond- 
ing dynamical patterns are drastically different. 

What conclusions apply to theories with less supersymmetry? In the 
simplest version of the Seiberg-Witten solution [1], Af = 2 supersymmetric 
QCD can be deformed by adding a mass term [i for the adjoint field. In 
the limit of large [i the theory flows to Af = 1 SQCD. At small fi the mass 
term for the adjoint field induces a Fayet-Iliopoulos F term in Af = 2 theory 
[36, 16], with £ proportional to fj, times some mass scale, such as A or quark 
mass. Thus, the deformation parameter \i translates, roughly speaking, into 
the FI parameter £. 

It is commonly believed that the behavior of supersymmetric QCD is 
smooth in /j,: the Abelian degrees of freedom of Af = 2 theory smoothly evolve 
into non-Abelian degrees of freedom of Af = 1 theory as we increase \fi\. 
While on the conceptual side our results provide an unambiguous evidence 
in favor of the smooth transition, dynamics-wise the emerging pictures on 
the opposite sides of domain lines separating domains I, II and III hardly 
look alike. In particular, light degrees of freedom are completely different. 

In addition we should note that there is at least one example of a non- 
supersymmetric model where the evolution is proven to be discontinuous, 7 
with a phase transition [5, 6]. And even in our basic Af = 2 model the 
crossover becomes a full-blown phase transition at N = oo. 

To conclude, we would like to comment on the recent paper [38]. In this 

7 On the other hand, analyses [2, 3] carried out in a nonsupersymmetric setting different 
from that treated here and in [5, 6] show no sign of the phase transition, while [37] exhibits 
a chiral phase transition on the way from Abelian to non-Abelian confinement. 
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paper it is argued that non-Abelian vacua with r > Nf/2 which support non- 
Abelian strings "dynamically Abelianize" in quantum theory. We disagree 
with this statement. As we demonstrated above, both the Abelian and non- 
Abelian regimes can be present in M = 2 QCD in quantum theory. They 
just occur in different domains of the parameter space and are separated by 
crossovers. 
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Appendix: 

Global symmetries of the CP (AT — 1) model 
with Z 2 AT-symmetric twisted masses 

First, let us outline the M = 2 CP(iV — 1) model with twisted masses [21] 
in one of a few possible formulations, the so-called gauge formulation [39]. 
This formulation is built on an iV-plet of complex scalar fields n l where 
i = 1, 2, N. We impose the constraint 

n\n l = \. (A.l) 

This leaves us with 2N — 1 real bosonic degrees of freedom. To eliminate one 
extra degree of freedom we impose a local U(l) invariance n l (x) — > e %a ^n % (x). 
To this end we introduce a gauge field which converts the partial derivative 
into the covariant one, 

d^V^d^-iA^. (A.2) 
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The field is auxiliary; it enters in the Lagrangian without derivatives. 
The kinetic term of the n fields is 



L = \ |V M nf 



(A.3) 



The superpartner to the field n l is an Applet of complex two-component 
spinor fields £ l , 

Ha- (aa) 

The auxiliary field has a complex scalar superpartner a and a two- 
component complex spinor superpartner A; both enter without derivatives. 
The full M = 2 symmetric Lagrangian is 



9l 



|vx| 2 +cVv^e+2^ 



a — 



V2 



\n 



i|2 



+ 



^ E (°- - 7=) 4^ + ^nj (X R C L - + H. 



c. 



(A.5) 



where rrii are twisted mass parameters. Equation (A.5) is valid in a special 
case when 



N 



i=l 



We will make a specific choice of the parameters rrii, namely, 



(A.6) 



(A.7) 



where m is a single common parameter. Then the constraint (A.6) is au- 
tomatically satisfied. Without loss of generality m can be assumed to be 
real and positive. The U(l) gauge symmetry is built in. This symmetry 
eliminates one bosonic degree of freedom, leaving us with 2A" — 2 dynamical 
bosonic degrees of freedom inherent to CP (A" — 1) model. 

Now let us discuss global symmetries of this model. In the absence of the 
twisted masses the model was SU(AT) symmetric. The twisted masses (A.7) 
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explicitly break this symmetry down to XJ(1) N x , 

n e -> e^n\ & - e**& £ - e^l , € = 1,2,..,JV, 
(7 -> fj , A RjL -> A fi)L . (A. 8) 

where are N constant phases different for different I. 

Next, there is a global vectorial U(l) symmetry which rotates all fermions 
^ in one and the same way, leaving the boson fields intact, 

t l R - e*dn &->e*&, £=1,2, ...,7V, 

n e ^ n e , a^a. (A.9) 

Finally, there is a discrete Z 2 at symmetry which is of most importance 
for our purposes. Indeed, let us start from the axial U(1)r transformation 
which would be a symmetry of the classical action at m = (it is anomalous, 
though, under quantum corrections), 

& - e*&, Si-e^, £=1,2,...,N, 
X R -> e^Afl, A L ^e"^A L , <r-*e 2i V, 

-> n e . (A. 10) 

With m switched on and the chiral anomaly included, this transformation is 
no longer the symmetry of the model. However, a discrete Z 2 n subgroup sur- 
vives both the inclusion of anomaly and m/0. This subgroup corresponds 
to 

7* = ^, k = l,2,...,N. (A.ll) 
with the simultaneous shift 

£^£-k. (A.12) 



In other words, 



Kr 




Xr - 




n e - 


+ n e ' k . 



(A.13) 
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This Z 2 n symmetry relies on the particular choice of masses given in (A. 7). 

The order parameters for the symmetry are as follows: (i) the set of 
the vacuum expectation values {(n 1 ), (n 2 ), ... (n N )} and (i) the bifermion 
condensate (£jj^£). Say, a nonvanishing value of (n 1 ) or {£ R e £ e L ) implies 
that the Z 2 n symmetry of the action is broken down to Z 2 . The first order 
parameter is more convenient for detection at large m while the second at 
small m. 

It is instructive to illustrate the above conclusions in a different formula- 
tion of the sigma model, namely, in the geometrical formulation (for simplic- 
ity we will consider CP(1); generalization to CP(A r — 1) is straightforward). 
In components the Lagrangian of the model is 

C C p(i) = G{d^dy-\m\ 2 <p^+ l -^[d R ^ L + ^ R d L ^ R ) 
. 1 - <t> ] <t> 



i 



X 

% 

X 
2 



(A.14) 



where 



and 



X = l + ^0, G=^— 2 (A.15) 

% X 2 

The Z 2 transformation corresponding to (A. 13) is 

^L^-^L- (A.17) 

The order parameter which can detect breaking/nonbreaking of the above 
symmetry is 

iM) Ski (A - 18 » 

Under the transformation (A.17) this order parameter changes sign. In fact, 
this is the central charge of the Af = 2 sigma model, including the anomaly 
[31]. 
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Now, what changes if instead of the M = 2 model we will consider non- 
super symmetric CP(N — 1) model with twisted masses? Then the part of 
the Lagrangian (A. 5) containing fermions must be dropped. The same must 
be done in the Z 2 order parameter. As was shown in [5, 6], now at m > A 
the Z 2 symmetry is broken, while at m < A unbroken. A phase transition 
takes place. 
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